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SYNOPSIS 

Tjie present analysis is aimed at studying the 
growth of a debonding crack in a composite. The model 
considered is a single fiber embedded in a cylindrical 
matrix with tensile load being applied to the fiber end. 
The fiber is assumed to have been bonded to the matrix 
everywhere except in a small region. Making use of the 
general solution for axisymmetric cylinders, a system of 
linear equations is developed. These are then solved for 
stress and displacement fields. The results have all 
been presented in a non-dimensional graphical form. 



1. INTRODUCTION & LITERATURE SURVEY 


1.1 Introduction 

The word composite means 'made up of two or more 
distinct parts'. From the earliest uses, the goals for 
composite development have been to achieve a combination of 
properties not achievable by any of the elemental materials 
acting alone; thus a solid could be prepared from constituents 
which, by themselves could not satisfy a particular design 
requirement. The struggle to develop new and better composites 
has not been limited to mankind, even nature has ventured 
into it and produced bamboo with excellent impact strength and 
buckling strength and bone with strength and lightness, to 
name only a few. 

Some examples of man-made composite materials are 
cermets, glass fiber reinforced plastics, reinforced concrete 
and plywood, cermets (composites of ceramic and metal) are 
used in cutting tools, turbine parts and as nuclear fuel 
elements and control rods. Glass-fiber reinforced plastics 
are being used in air-conditioners and humidifers, fire 
extinguishers, blood testing equipment. Steel fibers with 
matrix as unsaturated polyesters are used in chemical plants 
for polyester piping, ductwork, fans etc. and in building 
industry as thermosetting decorative sheet. And to top it 
all, composite materials have found their way in the more 



2 


strenuous and demanding roles like transportation . CFR 
epoxies have been used in DOT Fairchild experimental safety 
car/ while DC~lO and Lockheed L-1011 airplanes have used a 
hetrocyclic (aromatic) polymer named PR9-49, Looking at the 
tremendous studies made in the area of composite materials, 
the day does not seem far off when metals will be completely 
replaced by composites, possibly even in coinage. 

1.2 Literature Survey 

The majority of composites used in structural 
applications fall in the f iber-reinf orced category. The basic 
concept behind fiber reinforcement is the production of a 
two-phase composite structure in which deformation of the 
matrix is used to transfer stress, by means of shear tractions 
at the fiber matrix interface, to the embedded, high strength 
fibers. Provided the length of the fibers is sufficient, 
they should then be constrained to take up the same deformation 
as the matrix over the greater part of their length and thus 
effectively reinforce the matrix. In addition, the presence 
of the fibers retard the propagation of cracks and thus 
produce a material which is tough as well as of high 
strength. The present literature survey is concerned 
principally with the basic concepts of fiber reinforcement 
and the discussion of various theories proposed to predict 
the, thermoelastic properties of fibrous composites. All 
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these theories can be classified as belonging to either one 
of the following: netting analysis, mechanics of materials, 
self consistent model, variational, exact (within the context 
of classical elasticity), statistical, discrete element, semi- 
empirical methods and theories accounting for microstructure. 
Certain assumptions common to most of these theories are: the 
matrix and fiber are linearly elastic, homogeneous and free 
of voids; there is complete bonding at the interface of the 
constituents and there is no transitional region between them, 
the ply is initially in a stress free state and the fibers are 
regularly spaced and aligned. 

Netting analysis assumes that the filaments provide 
all the longitudinal stiffness and the matrix provides the 
shear and transverse stiffness and the Poisson effect. This 
method provides acceptable values for longitudinal properties 
and low values for transverse and shear strength. In mechanics 
of materials approach, the ply's averaged stress-strain 
temperature state is expressed in terms of the stress-strain 
temperature states of its, constituents using displacement, 
continuity and force equilibrium conditions. Substitution 
of the ply averaged stress -strain -thermal state in the 
appropriate definitions for thermoelastic properties yields 
the desired results. 

HilltkQmodeled the composite as a single fiber 
embedded in an unbounded homogeneous medium which is 



4 


macros copi cal ly indistinguishable from the composite. The 
medium is subjected to a uniform loading at infinity which 
induces a uniform strain field in the filament. This strain 
field is used to estimate elastic constants. Hill's method, 
coming under the self-consistent model approach gives reliable 
values at low filament volume ratios, reasonable values at 

intermediate ratios and unreliable values at high ones. 

t 

Rosen has used variational method to obtain bounds 

A 

in the ply's elastic properties. The minimum complementary 
energy theorem and the minimum potential energy theorem have 
been used to obtain lower and upper bounds respectively. 

An exact method of analyzing filamentary composites 
consists of assuming that the fibers are arranged in a regular 
periodic array. The resulting elasticity problem is solved 
either by a series development or by some numerical scheme. 
Sarma^has solved for stresses in a f iber-pullout model using 
exact method. Statistical methods assume that the composite 
is a homogeneous solid whose thermoelastic properties vary 
randomly with position. One fluctuating term is added to the 
field equations. The resulting equations are statistically 
averaged to give the desired results. 

Discrete element methods use a finite element 
scheme (usually a constant-stress triangular element is 
employed) to solve for elastic constants. Semi -empirical 
methods tend to retain the physical complexities of the 
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problem while relaxing the mathematical rigour. Bolotin has 
proposed a theory accounting for microstructure assuming the 
fiber to be more deformation resistant than the matrix; the 
distance between fiber and matrix and thickness of fiber being 
small compared with the characteristic dimensions of the body 
and comparable with the distances over which the functions 
characterizing the stress -strain state of the reinforcing 
elements change appreciably. This equations resemble the 
ones for a Voigt-Cosserat medium, that is, a medium for which 
the stress tensor is not symmetric. 

Research work on composite materials has primarily 
been centered around models assuming no solubility or no 
reaction at the interface. Also properties were not associated 
with the interface, per se £4] For example, the term 'inter- 
facial shear strength' was often applied to the stress in the 
matrix immediately adjacent to the filament. A further 
assumption was made that the interface was stronger than the 
matrix so that the matrix flow limited load transfer from and 
to the filaments. More practical systems, now in vogue have 
used higher strength matrices, causing the failure path to 
move to the interface and this has led to the development of 

theories for weak interfaces. 

L23 

Cooper and Kelly^have divided the mechanical 
properties of composites into those affected by the tensile 
strength of the interface cr , and those dependent on the shear 
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strength For longitudinal tensile loading, they conclude 

that the tensile strength o± the interface is not critical, 
but the interface shear strength controls the following 
properties, critical or load transfer length i.e. the filament 
length required for the longitudinal stress in the filament 
to reach its fracture stress; composite fracture under condi- 
tions of fiber pullout and the deformation of the matrix in 

C3] 

fracture. Ebert and Gadd have modeled the composite by 

A. 

considering a typical fiber-matrix unit to have a cylindrical 
configuration with a core of matrix and an outer sleeve of 
fiber. They conclude that the transverse stresses generated 
in the comoosite reach their maximum intensity at the fiber 

ey 

matrix interface. Lawrence has discussed the distribution of 

A 

load and shear stress along the fiber length and the condi- 
tions for fiber pullout. He concludes that the maximum fiber 
load depends on length of embedded fiber and ratio between 
the shear strength and the frictional shear strength of the 
interface and that the development of debonding of the fiber 
from the matrix can have a marked effect on the maximum shear 
stress developed at the interface. 

1.3 Sc ope of the Present Work 

The present work has grown out of the realization 

( 

that the majority of work done in the past has started 
either from the assumption of a perfect or no bonding at 
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the interface. The most important cases of limited strength 
bonds, or regions of good bonding mixed with regions which 
are unbonded have received little attention. It is the 
problem of regions of good bonding mixed with unbonded 
regions that interests us here. This problem is an extension 
of results obtained from Sarma'^ljThe elastic solutions as 

tsl 

obtained by Murthy has been used for numerically evaluating 
the stress and displacement fields. 
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2 * ELASTIC ANALYSIS 

2.1 Choi ce of the Model 

In the present problem,, the model considered is a 
single filament embedded in a cylindrical matrix. The choice 
of the model was primarily dictated by its resemblance (in 
geometry) to actual composites . Both the fiber and the matrix 
are assumed to be linearly elastic. The matrix is clamped 
on the lateral surface, clamping is also done on the top 
surface of the model (see Figure 2.1). The load is applied 
to the bottom end of the fiber. The debonding is assumed to 
have occurred over a distance c at a place XDB units away 
from the top end. 

Lame has solved the problem of determining the 
stress and strain fields in a circular elastic cylinder under 
various boundary and loading conditions. His solutions, 
however, assume a plane strain condition. For cylinders 
whose length is comparable to diameter, the plane strain 
condition is no longer valid . 

2.2 Analysis of E lastic cylind ers 

The stress displacement equations for an elastic, 
homogeneous medium in the presence of axisymmetric loading 
can be written in the form? 
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O' = 

cr = 0 

(1) 

r0 

ez 

a =s 

r 

(X + 2,) x<g + ^) 

(2) 

a e 

<> + 2,)H + A(|| + ^) 

(3) 

a = 

Z 

<*+ ’•*>£+ ><£ + ?> 

(4) 

a = 

rz 

•* ^ + n> 

(5) 


where X and \x are the Lame * s constants . 

The equilibrium equations to be satisfied by an 
axisymmetric problem are: 


^ g rr 3> °rz 
~b r D z 
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( 6 ) 


e c rz 

"a r 




__ 
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°rz 
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★ 

+ Z 


o 


(7) 


ie 'k 

where R and 2 represent the body forces per unit volume, 
in the r and z directions respectively. 

Substitution for cr , cr and cr from Eqs. (1) to (5) 

X. 2TZ 

in Eqs. (6) and (7) yields: 

( > + 2 n ) + . 2_, + p. + ( A + n ).-|^u + r* = o 

or r " 

( 8 ) 


( /.+ |i ) ~ ^ + 2 ) 


2 

o w a jx ^ w 


?f 


w 


jr^z Tjz- ' ' - - ' " r Dr + ^5 r 2 


+ z = o 


( 9 ) 
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Elimination of w from Eqs. (8) and (9) leads to the 
fourth order governing partial differential equation in r and 
z as; 


3 1 "i 

r 5"r 


r + 




2 

r) u 


+ 


o * 

JR 

‘ 2 
2) z 


+ 


( \ + 2\x) 3r r }r 


r R 


( ' *t ) 

P- ( + 2 P-) 


-v 2 * 

•J Z 
?} z - ? } r 


0 


( 10 ) 


For the cases involving constant body forces, Eq. 
(10) reduces to; 




r 


•}r 



u 


0 


( 11 ) 


Thus once the solution of Eq. (11) is known, w can 
be calculated using Eq. (8). The calculation of stresses and 
strains is then a simple exercise. This has already been 
done by Murthy. His results form the basis for the following 
analysis . 


2.3 Displacement Equations for Fiber 

f f 

The displacement equations for u and w (radial 
and axial displacements for fiber) are taken as; 

f , 3 2 3 3 „ 

u = (S^r + S 3 r + S 5 rz + T 1 r z + T 5 rz ) 

MB § 5 

+ XL r J (6 r) [s, cosh 6 z + T, sinh 6 z] 

« o m 1m m lm m 

m=l 



11 


NSUB cosh p z R sinh {3. z g 
^ ^^rAs v+r-rATfj Jjp.r) 
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-cosh p,l lk sinh (3 1 lk J 1 H k 
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+ [rj, (6 r) 
i m 


'•f + 2 ^f 2J o ( m r) l: 


iS, sinh f -. z 
. 1m - m 


£ NSUB r sinh j ,z jl cosh , z g p 

+ x lm costl m Z ' ~ ; cosin j-'^l S lk + sinh"T^l T lk j J o^‘ k r ^ 

(13) 

where ip s and • v 's are eigen values obtained from J 1 ( r a) = O, 

JL k 


J o ( - j ) -° - a y n -T- 


Knowing u and w , the stresses can be calculated 
using equations (1) to (5). 


2.4 Displacemen t Equat ions for Matr ix 

The displacement equations for u m and w m (radial and 
axial displacements for matrix) are taken as 

u m = \ E-j r 3 + E-r + E c -rz 2 + F n rz 3 ; + ^ rY i'-.r) >E^. cosh 5\z 

t. 1 3 3 1 ) j=l ° J 1 3 J J 


£ , NSUB . cosh p z & sinh & z 

4 - F_ . sinh C, .zv + ^ _ _ * e" + 


1V+ rirfh-r-T t-J-, (fi r) 


cosh 0 1 Ik sinh hj, 1 lk \ l u *k 
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r W w> 
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+ ^ — j— — t — F 1 + — ~ — T-T — r— F _ f sin x z 1 

(W ln VV 3n J Xn 1 


( 14 ) 
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m . 


The corresponding expressionf for w is: 
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Knowing u™ and w m , the stresses can be calculated 
using equations ( 1 ) to ( 5 ) . 

2.5 Boundary Conditions (Exact) 

F or Matrix ; 



u(l, z) 

0 

(16) 


w(l, z) = 

0 

(17) 


*\<r* = 

0 

(18) 


C" ( r , o) = 

0 

(19) 


^ rz (r ' X) 

= 0 

(20) 


- r rz (r ' o) 

= 0 

(21) 

For 

Fiber; 




<- 2 (r, o) 

= 0 

(22) 


cT* ( r , 1) 

= £ 1 

(23) 


(r, o) 
rz ' 

= 0 

(24) 


cr (r, 1) 
rz 

= 0 

(25) 

At_ 

the Interface; 



(i) 

For bonded region 

(z < XDB or z (XDB + c)) 

/ 



u (a, z) 

= u (a, z) 

(26) 


m / s 

w (a, z) 

= w (a, z) 

(27) 


^ 2 (a, ,) 

= rr f (a, z) 
rz 

(28) 


z) 

f f V 

= cr (a, z) 

(29) 
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(ii) For debonded region (XDB z < (XDB + c)) 


( a / 2 ) = 0 


cr' r (a / z) = 0 


<r* z (a. z) = o 


cT™ (a, z) = O 
rz 


Condition for Force Equilibrium 


z=l 

f c — (a, z) dz 

rz 

z=o 


2 qa 


(30) 

(31) 

(32) 

(33) 


(34) 


2.6 Determination of constants in Equations (12) - (15) 

The constants in the equations for u^, u m , w m 

etc. are determined by the application of various stress and 

displacement boundary conditions. Truly speaking, the 

enforcement of proper boundary conditions would have required 

f m 

all summations in u , u etc. to have extended to infinity and 
this would have entailed the solution of an infinite matrix 
for constants. However, an approximate solution with a finite 
matrix can be obtained by truncating the infinite series. 

Even with this simplified solution, one cannot 
satisfy the boundary conditions at the interface exactly 
due to the change in the conditions between bonded and 
unbonded regions. A collocation technique was therefore 
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employed to satisfy the boundary conditions in the debonded 
region. For convenience, these points are chosen at uniform 
intervals. Thus the approximate boundary conditions for 
debonded region become; 

^(a,XDB+^~c) = 0 (35) 

cr f ( a, XDB + ■ c) = 0 (36) 

r k-1 


c rz a ' 

XDB 

M-l 
+ k-1 

c) 

= 0 

(37) 

cr f (a, 
rz 7 

XDB 

M-l 
+ k-1 

c) 

= 0 

(38) 


where XDB and c represent the location and extent of debonded 
region (Figure 2.1), k represents the number of subdivisions 
within the debonded region and M is an integer variable such 
that 1 M £• k . 

Obtaining the set of equations from the above 

boundary condition is fairly straightforward. What one does 

£ m 

is to substitute the expressions for u , u etc. in the 

various boundary conditions. This leads to an expression 

either in z (for equations on z like u m (l, z) = O) or in r 

£ 

(for equations on r like dr (r, 1) = f . ) or in terms of constants. 
This expression in z (or r) is now broken down into simpler 
expressions and each of them is expressed as a half -range 
Fourier Sine series (Fourier-Bessel series for expressions 
on r). Regrouping the expression, one obtains a Fourier=Sine 
(or Fourier-Bessel) series on the left hand side. One can 
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follow a similar procedure on the right hand side. A term by- 
term comparison of the series on the two sides leads to the 
desired set of equations. For example, consider the boundary 
condition u m (l, 2 ) =0. Substitution of the expression for 
u m in the above boundary condition gives: 


u m (l, 2) = (E 1 + E 3 + E 5 2 2 + F 1 z 3 ) 


NF 


+ ^ 


j=l 


rY ( 6 .) Se” . cosh t'S.z + F. . sinh S .z 
0 n 3] 3 3 \ 
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E ik + iiHF'fjT F ik ^ J i (i k } 


MF 
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+ tSi rJ (6>. ) -? E- . cosh 6 .z + f 5 . sinh <S .z r 
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O J <1 3-J J 1 3 


7 
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NF 


C r ^ K, V 

+ ^ E n „ + E 0 _ '1 COS <^Z 4- (F- + l sin '\2 


n=ir in 3n J 


n 


In T 3n 


n 


O 


Replacing each function of z by its Fourier equivalent, 
we obtain after the regrouping of terms: 

NF - 


cr.. 

> 




I (E, + E-)xr + E.X 3 + F X~; + s. rY ( t. ) )E'Xf f;.X°. < 
Ll 3 m b m 1 m o J ) 3] r.ii 3j mj - 


+ \ + E ^3 X 1 + [ F fn + F 3n>n 

n=i 1 . ^ 


NSUB <^cosh z sinh (5 z k 7 

7 ^ j coih~^~T E lk + sinh~^T F lk sin ^ 


0 
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Thus equating each term of the Fourier Sine series 
to zero, we gets 


NF 


1 3 ^ XMJ X n cS 8 "J 

(E, + E_)X + EJC + F.X + Hz. rY (,h.)>E_.X . +F..X ( 
1 3m 5m 1 m ^ o j ' 3j raj 3j mjj 
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X^l ! cos h ’^>^1 sinh p>^.l Ik I 1 
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Similarly equations can be obtained for the 
remaining boundary conditions. For clarity sake, the various 
terms have been grouped under different categories. They are: 
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5. 1 5 3j 3j 

c h 

s ij' T ij. 


= 1 


Si / S 0 , Sr , T 1 / Tr, S 01 ..S. 


< 


A _§ 


cf 


x l # i 5 ' w 31 # * w 3m' T 31** T 3m' S ll** S lm' 


r 

c 


T ii-- T i m s 


<Si , 

c 1 


, S 3 , S 5 , Tl , T 5 , <S° j ) T ' ( (T3 1 ) T , (S“.) T , 4^ 


6 S T 

3j 




4 ,T;T 

ly ^ 


The other vectors are: 




,o \T , o .t 


(EKNWN) = y|E 1 , E 3 , E 5 , F 1# (E^) , (F 3j ) x , (E^)\ (F~ . ) 1 , 


(4M f +l) 


T 1 X T 1 T 1 _ 

(E' ), (F"l . ) T "? T c 

3j 3j Jt / 
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(SALPH) = 
( 2N f ) 


' (s"' S ‘ ) F (t >v )^' F 

/ iS lj ; 7 U lj ; y 


(EALPH) = ; (E\ . ) , (Fi . ) 


( 2N f ) 


(SBETA) 

(2k) 


T , T 


i (s ij )T ' < T i j )T ; T 


(EBETA) 

(2k) 


■; (E ij )T ' (F ij )T - T 


(SCONST) 
( 4M£sf-l ) 


Load vector associated with fibers 

T IT T T 1 nn T 

; (0) , (PHI*G m ) , (0) , (0), (Pl*G m )' 


Thus the equation u ( 1 , z ) =0 can be written in the 
matrix form as; 

lx 1 X 1 X 3 X 2 X 6 m ( 0 ) T ( 0 ) T x 7 X 8 \ * 

i m m m m mn n v ' mn mn J 

> E 1' E 3' E 5 ' V (E 3j> T ' <4 )T ' < E lj )T ' < E ij> T ' < E 3j )T ' < E 3j )T,T 


J, ( P . )X . J, ( p . )x . . 

1 ' J i ' -i l ' -1 ' -ro -i i 


, j L_J iiUL) (__± L — HD.) 1 \ ( S f ) T CT I ) T ! - o 

+ L cosh .1 ; 'sinhft.l ( ' ' S 1 j ' ' ^ T lj ; 1 ° 

V 3 f J 


(B ) (EKNWN) + (B ) (EALPH) + ( B ) (EBETA) = (0) 

M M M 


where (B ' ) = \ X* X* X 3 X 2 X*L 
M , m m m m mn 


m o T o T x 7 x 8 ’! 

n mn mn -, 4 
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it ft <■ 

(b m ) = (x® ^ ), 

M mn mn 

9 10 

... X .e. X . e-, 

(B M cosh f- .1 sinhj-3 .1 

t j rj 


where 



The equations to be satisfied on the fiber, matrix 
and the interface are written in a form similar to the above as 
shown in the Appendix II. 


Solution of Equations ; 

The equations can be written in the following form. 

(i) Equations on fiber; 

(aJ) (SKKTWN) = (A^) (SALPH) + (aJ) (SBETA) + (SCONST) 


l -^(4M £ +1 )x(4m £ +1 ) 


(4M f +l)x2N f 


(4M £ +1 )x2k 


(4M f +l ) 


(i) 


(ii) Equations on matrix; 


<B*) (EKNWN ) = (B^) (EALPH) + (B^) (EBETA) 


(4M £ +2N £ )x(4M f +2N £ ) (4M f +2N f )x2N f 


(4M £ +2N £ )2k 


(ii) 


(iii) Equation on interface (Bonded region); 

(A^) (SKNWN) + (a|) (SALPH) + (a|) (SBETA) 

4N £ x(4m £ + 1) 4N f x2N f 4N f x2k 

= (B*) (EKNWN) + (B*) (EALPH) + (3*) (EBETA) 

rtM — f 4N £ x(4M £ +2N £ ) 4N f x2N f 4N f x2k 

(iii ) 
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(iv) Equations on interface (Debonded region): 

(aJ) (SKNWN) + (a|) (SALPH) + (A?) (SBETA) 

r *2kx(4M £ +l) (2kx2N f ) (2kx2k) 


(B 7 ) (EKNWN) + (B®) (EALPH) + (B S ) (EBETA) 

M MM 


= < 0 ) 


= o 


'2kx(4M £ +2N £ ) 


2kx2N, 


2kx2k 


From equations (i) and (ii) 


(iv) 


(v) 


(SKNWN) 


= (k.^)" 1 (A^) (SALPH) + (A*)" 1 (A^) (SBETA) 


(EKNWN ) = 


1 -1 

+ (Ap) (SCONST) 


(B^)" 1 (B^) (EALPH) + (B^r 1 (B®) (EBETA) 


(vi) 

(vii ) 


Substituting the values in (iii), (iv) and (v), one 
gets after rearrangement: 


(AF412 

5. 

+ V 

4N £ x2N £ 


(AF712 

+A ! } 


■(BM412 (AF413 


(O) 


2kx2N, 


5, 

+B) 

M 

4N £ x2N £ 


(O) 


1 2kx2N f 2kx2N f 


(BM712 

+B®) 

M 

2kx2N J - 


6 , 

+Ap) 

4N £ x2k 


(AF713 
9 N 

+Ap) 

2kx2k 

(O) 


2kx2k 


(4N J= +4k)x(4N,+4k) 
r r 


’ ( BM4 1 3 1 f (SALPH) 1 

««£ 1 1 

4N £ x2k 


(O) 


2kx2k 


(BM713 

+B^) 

M 

2kx2k 


(EALPH) j 

t 
\ 

S ) 

I (SBETA) | 


(EBETA) 


J 


J 


(4N f +4k) 


-(AF41)x . 

t 

(SCONST) j 

(4N £ ) j 
! 


-(AF71)x 

i 

(SCONST) 

(2k) 

(O) 


(2k) 




(4N jC +4k) 

1. 1 L *' .(vi 4 ilil|t 

CEH Ik 1 -l ' 'yRARV 




59493 
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(AF412) 

= 


<4 rl 


( BM4 1 2 ) 

= 

<b m> 

(B M rl 

(b m> 

(AF413) 

= 

( 4 ) 



(BM413) 

= 

(B «> 

(B M rl 

(b m> 

( AF 712) 

= 


( 4 } ~ 1 

(a f 2) 

(BM712) 

S= 

< b m> 


(b m> 

(AF713) 

= 


(Ap ) 

< 4 > 

(BM713) 

=S 

(b m> 

( bV 1 

M 

(b m> 

(AF41) 

= 

<4> 

(Ap )" 1 


(AF71) 

= 

<4 ) 

( 4)- 1 



The equation (viii) is a set of (4N f +4k) equations 
in (4N^+4k) unknowns. Thus it can be solved for (SALPH), 

(EALPH), (SBETA) and (EBETA) . Then (SKNWN) and (EKNWW) can 
be determined using (vi) and (vii). Knowing all these constants, 

£ m 

the values of u , u , etc. can be found by simple substitution 


in Eqs . (12) to (15) . 
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3 * DISCUSSION & CONCLUSIONS 

3.1 Numerical Computations 

The following data are used to get the solutions 
for a hypothetical case: 

Poisson’s Ratio for Fiber = 0.25 

Poisson's Ratio for Matrix = 0.27 

Ratio of Elastic Moduli! of Fiber to Matrix = 20 

Non-dimensional Length = 2.0 

Non-dimensional Radius of case (Matrix) = 1.0 
Non-dimensional Radius of Case (Fiber) = 0.4 

Applied Load (Non-dimensionalized ) 

= Ratio of Load Applied at the Fiber End to the 
Shear Modulus of Matrix 

= 0.00001 

Non-dimensional Location of Debonded Region (XDB) = 1.2 

Non-dimensional Extent of Debonded Region (C) =0.5 

For convenience, the computations were carried out 
in two stages. In the first stage, the coefficients of 
Fourier-Bessel series were evaluated using exact formulae 
wherever they were available and resorting to numerical 
integration (Simpson’s rule with 100 interval) elsewhere. 

The results so obtained were used as input for further 
computations . 
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In the second stage, the infinite set of equations 
(Appendix II) are solved by considering only three terms in 
the solution. The LU-decompcsition method is then employed 
to obtain the variables of the equations . The displacement 
field is then obtained by using equations (12) to (15). 
Stresses can be computed by using equations (2) to (5). 
Appendix III gives the flow chart for computational work which 
has been used to develop the programme for DEC 1090. 

3.2 Results 

Figure 3.1 shows the plot of axial stress with 

length. It is observed that k attains a maximum at the 

z 

middle of the length and goes to zero at the two ends. 

Figure 3.2 depicts the variation of shear stress 

with length. The shear stress CT is zero at z = 0 and 

rz 

reaches a maximum value at z = 1/4. it tends to zero at 
z = 1/2 and remains constant beyond this value. 

The variation of radial displacement with length is 
shown in Figure 3.3. It is observed that the displacement 
increases continuously till z =-1/2. It then keeps falling 
down till the end of the debonded region and then remains 
nearly constant. 

Variation of the axial displacement with length is 
shown in Figure 3.4. It is observed that the axial dis- 
placement increases monotonically with increasing z. 
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Figure 3.5 is drawn to show the variation of /T' 
with length. The radial stress reaches a maximum in the 
bonded region and becomes zero as expected in the unbonded 
zone . 

3.3 Conc lusions 

Following conclusions can be drawn on the basis of 
this limited study. 

1 . The axial stresses are of a much higher order 

of magnitude than the shear and radial stresses, which implies 
that the load transfer is mainly through axial stresses. 

2. From Figures 3.2 and 3.5 it can be seen that 

the stresses and *t* reach their maximum values in the 

rz r 

upper region of the length (0< z <^1/2) and hence there is a 
liklihood of crack initiating in this region , rather than in 
the lower portion . 

3.4 Suggestions for Future Work 

1. It is suggested that this analysis be carried 
out by taking a few more terms in the infinite series to study 
the convergence of the solution and also give a better picture 
regarding the variation of stresses and displacement along z. 

2. The model can be improved by considering multiple 
fibers in a given matrix which may result in redistribution 

of stresses in the debonded zone. 
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APPENDIX It EXPANSIONS 


Expansion 


Remarks 


1 — G J ( r ) 
v \\ mom 

r = T g 2 J r) 
mom 


r 2 = V g 3 J ('■ r) 
mom 


i 


3 v~- 4 

r = / G J (:• r) 
m o m 


1 ~ ( r) = g j (; r ) 
° n ; s mn o m ; 


G = 
m 

G 2 = 
m 

G 3 = 
m 

G 4 = 
m 




2 /J l (h m } J r 2 j o ( 'n, r)ar 

2/J l (< m> J 1 “o' m r Mr 
o 


1 1 ( 

\ r) = 
n 

= G 6 J ( 

, mn o 

*' r) 
m 

G 6 = 
rnn 

H 

0 

( -- r) 
n 

= ~ G 7 J 
mn o 

C r) 
m 

G 7 = 
mn 

rip^r) 

= ; r g 8 j 
rnn o 

( ? r ) 
m 

G 8 = 
mn 

I 1 (- 

Si r )/ r 

= ( '~ G W J , 

_ ( r ) 
o m 

G 9 = 
mn 

J (} 

^ r) = 

G 10 J (l 

; r) 

G 10 = 

O 1 

n 

.,. ;i mn o v 

m 

mn 


r) - 
n 

21 G 1X J ( 

-) mn o 

, r) 
m 

G 11 = 
mn 

J 1 ( F 

' r)/r 
n 

12 

= G Z J (V. r) 

. , mn o m 

G 12 = 
mn 


2 /J l ( “m ) j rI l (= ' n r)J o ( i, r)dr 


?(£ ) ! r 2 } 

1 m 


2 **" 2 

2 /J l ( m ) -- r I i (>i n r ) J ^(^r)dr 


o m 


o n o' m 


o m 


= 2/J l C m> i 


2J (/-" )$ /( 2 =?- 2 )j /v } 
o < n' n/ Vr n m jvJ l^ m J 


2/j l ( m ) J rJ l ( fn r)J o ( '~m r)dr 
2 c 

2 /j l ( 'm ) j J l (r h r)J o (c ' m r)dr 
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rJ_(£_r) = / g^j (,s r ) 


o n 


,13. 

mn L o m" 


; 13 
y m n 


rot r) = ‘ G 14 J (<■' r ) g 14 

1 n ; mnV c m r; G mn 

G 15 
mn 

G 16 
mn 


= 2/J l ( ” m> ^ r2j o (f n r)J o (? m r)dr 
= 2/J l ( V -< r2 J 1 <i n r)J 0 (? m r)dr 


J-, O' r) = g 15 j (f r) 
1 n -.1 mn o' m ' 


rY o ( ^n r) “ " G 16 J (<, r) 

on - - mn o' m ' 




r)dr 


m- g ““I' n o- m 


2/J l ( ‘V | 


r Y,(J r) = ?. o 17 J (t r ) G 17 
x n — mn o m mn 


2/J l ( "m ) j r 2 YiC_ r )Jo(? r)d r 


n o v m 


Y 1<‘V> 

= G 18 J (c r) 

- mn o m 

G 18 

mn 

= 2/jJ 

1 

a m ) J rY 





o 

V'n r > 

= , G (c r) 

- mn o K n J 

G 19 

mn 

= 2/j^ 

1 

<: 'm ) 2 rY < 





o c 

J o<V> 

= / x n j ( r) 

/_ m o' m 

, n = 

~ v m 

l 1 

if n = m 


' V- i 

1 o 

otherwise 

i = ; 

1 • / 

X„ sin < z 
n n 

yl = 
n 

2(1 - 

( -1 ) n )/'vrn 


-j V JL /U ^ 

in o 


m 


sin ^ z 


n 


n 


X n = =21/ (n T i ) ( =1 ) n 


2 3 . 

z = -■ x sm z 
-- n n 


= “ 21 2 / (n-fr) - 21 2 /nV 3 (l-(-l) n ) 


3 -- 4 

z = /- X sin n z 

u~i n n 


4 5 

2 = X„ sin z 

>s . ; n n 


X n = -21 3 (-l)”/m; + S/nV(.l] n 
x 5 = 


n 


”21 4 [(“l) n / nT T + 12/n 3 T. 3 
+ 24/n 5 TT 5 (l-(-l) n )/ 
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COS .< . Z 

J 


X sin 2 
n J n 


, 5 0 if (n-j) is even 

y^n . =: 

4n/(n -j ) sT if (n=*j) is odd 


cosh h .z = 


.7 


.. X. sin x 2 x 7 . 
n J n nj 


2n~' 


(nV + <1+< - 1>n cosh fjl) 


sinh i .z 

J 


8 Q 

x ni sin -z X 8 . 

n J n n j 


2n i' n 

7TT — TT (- 1 ) 

(n -rr 2 + f 2 i 2 ) 


cosh ( -jZ) = 


,,9 . q 

X • Sln x • 
n J n nj 


2n * * / n « 

~ 2 2 ~2~2~ (f+C-l ) cosh j?- .1) 

(n rr + - r j 

! J 


sinh(j'.z) = 


X^? sin x 2 x 1 ? 
< n J n nj 


2n_El_ , -j vn 

,2 2 , 2 , 2 . ( “ 1) 

(n it + p .1 ) 
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APPENDIX II; DETAILS OF EQUATIONS APPLYING BOUNDARY CONDITIONS 


( i ) E quations on Fiber 

< “4(3. \ r +4) ~ ^ «=( t +2) 

l/F'-npi r-°«> V 2 Vi> Vr7pn- G » ) (0) <0> 


(3/ : -+4) G . G . (3,'j-+4) 

2 + rf^rr 1 <°> 2 , p^‘r^vrr-^y 


14' m -? i 

+ G mj } (0) 1 (SKNWN) = j (2 f v F -^(-f~f^ 7 y) Co) » (SALPH) 


- [ (2 • F .'j ( c osl l 1 '-n: )) (0 >; (SBETA > 


(A. 22) 


8( a r+2) o 1 4(3,\j.+4) q i ( *\ jr+2) -j 

? ( 7 - 3AG - -7 G ) '' (2‘A.G ) < (- 7 i-rr G j 

! - F ( \ f +l) m ( a £ + 1) m F fra F ( \ £ +l) m 

2), 91 4(3-44) 8(2 +2) 

- — - ~=rr 1 G ) , (- T _-f T — 1G" 3 + -rr-. l^G 

A f ^ ^ ( A £ +1 ) in 3 ( ^ ^ -hi ) m 

3(; ,+2) - 2 a 3 . ( 3 • +4) G 3 

7— -~4p~- 1G 3 - T -~ T 1 3 G X ) (-2 « >.(-1 ) J (t — — y v 

•'F ( > f +l) m ( 2 f +l) m F j < /- f +l) j I Q (<-<.) 


(3 2 f +4) 


o 'j 


hr) (0) 2 f cosh j 1( - TTT+iT" “h + G mj> 

J _ r j 


(3 • £ +4) /% ‘ 14 ■ 

2<„$..sinh .!(- 7 tux - 1 + G.) ! (SKNWN) = 


- F ""j 


j v ( ++D 


G 5 

r _2' lF - j ( -D 3 <i^ T ) (0)'! (SALPH) + 12 : F !ys“) 

Y . ( G^ G ) 'j ( SBETA ) + R PHI X G*)] 

3 3 ' 1 


2 


(A. 23) 
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(3, v +4) 3( +2) . 

O g!) (0) 


(0) (O) (o) ! ("Ji “ — 'j ) \ Tj — j? -s—c- '~ T 

F " ■ s f +1 ^ ra F v 4( ,\ £ +l) m 


( \ £ + 2 ) G_ 


G 7 . 


2 : tV( Tpr? Tp^fcjT + i o T 1 TT ) (0) 2 --F'j ( Trf+n 


( j-+ 2 ) G^ no 

-^1 + g 13 
- j 


G 6 . 


(SKNWN ) = • (O) 2 VF ^(- YT L 7) ) 1 < SALPH > + j’(°) 

1 j *- 


G 


11 


2 ' F ij ( ” : (SBETA) 

>' 3 


(A. 24) 


8( J f +2) ,. 2 , 


(3 A^+ 4 ) 


("7 TT7 1G^) (0) ^(T~rr 1G!) < ( 7 7 TT - G 


*F ' ( .•. £ +T) 


4 


m 


F v ( > f +ry '-FM ,. £ +l) m 


4( 'f +2) ,2 2 , . 3( 'f +2) 4 3A f 2 „ 2 , - , 

( ,\,+l) 1 G m ' F 4 ( ' +1) G m + ( a .+1 ) 1 G m ) (0) 


. ( - ,+2) G 6 . G 7 . 

2 ' F -V-) J <T-^ry:^hjy + ^-T> 


( ; +2 ) G 15 , - 

f mj „ 1 3 


2-'. .. .sinh . 1 ( *7 — 

’F J J ( - f +l) .. j 


( \ +2) G 1 ^ 

+ G ,) 2 i A. COSh. . l(“p— tr-r- — ~r * ^- 

- F 3 ‘ 3 ( .'.£+!) V j 


+ G*f) ! (SKNWN) = ! (O) 2 f X. (-l)-J(-. - m 

m 3 , . -F 3 I l v> ’j 


| (SALPH) 


G 11 . 

+ ! 2 fitanh ! : _.1 (~g"") 2 ; p ^(- ^—p-j—j) ] (SBETA) 


•F : j 


, 11 . 

* 3 ' mj 


^ 3 


(A. 25) 
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r . 8 (/i f + 1 ) 2.2 -• 

! 7- r r- T i-s rr^) (0) i 2 a 


:-f ' ( +1 ) T 


F v ( r Vf +i) 


<L^ -I 3 
F l ( ,-,,+1) 1 a 


4 ( a _+2) ,3 

t al 



2 


2 ' f _ cosh , . 1 ( 
F i 


( y<^+2) ( <A) a ^„("V,- a ) 


u. 


( A f +D 

^ j 1 0 ^ ^ 

( \ £+2 ) 


( A f + 1 ) 

' i 

( :\ f +2) 

J 1 C i a) 


o ( ' K j 


) (0) 


+ aj ( i..a) ) 
o J 


2 ( sinh . 1(7 t-v 

F 3 ( \ f +l) . • ' o ' j 


r ■; F 1 ^ ^ ) , 

j 2 ^(l-(.l) J )( r | . ^ ) (O) ! (SALPH) 


1-cosh y .1 

+ 2 f J- ( f .a) ( r-T-J-) 

* ’ F l i 1 cosh , . 1 


+ . Pl*G 


1 1 
'm I 


2 ; J- (o a) (-1) ) | (SBETA) 

(A. 34) 


(ii) Equations on Matrix 


K> < X m> <*£> < X m> <X mj> <°> <°> 


(Y Q ( > j ) )X®_. 1 (EKNWN) 


(~X^\) (- &“> ; (EALPH) + 


^1 0^ ) q ( i -j ) X0 ‘ 

( t J v ) („ X . ) ( — r~ — ^ — r ) ( “X .) i (EBETA) 

cosh m j smh mj . 


(A. 16) 
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f( m - 2 
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(0> 


(EBETA) 


(3 >\ +4) G 19 

T—~rr- + g 7 ) 

( > m +" j ">J 


— -1-y) (O) ~| (EALPH) + 


(A. 19) 
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f 8( > m +2) 2 1 4(3 ^ m +4) 3 


< A m + 2 > 


m ( 


G) (2 Vi G ) (- -r-iH-T-r- G' 

m m' +37 r 


2-^ m 2 1 4 ( 3 ,\ +4 ) _ 8(^+2) „ , 

1U "j x \ / m m -2JL V 

*/ ^ ,fT ' ' '** / i \ lG -f 0 v 1 G ) 

( >m +1} m ( ^ m +1) m 3 ^ ; 'm +1 ^ m 

-; ( 3 / +4 ) G G 

( - 2 V- 1)J wn^:-niT + i7b> <°> 


•m j o ' \j 


o j 


(3 A +4) /N m 

(2 ; ". cosh K 1 ) x ( •=■ j — - — S..— -.-l-J- + gJJ) (2o, sinh ,-.l)x 
J J ' /' m +1 >j m J J J 


(3 A m +4 ) 

-r— ' 4- C -r 


- L + G .) (2b. cosh ^.l)x(. 

< A + 1 ) j j 


(3A +4) G 1 ^ n _ 

' m mi 17. 

( .4 +D c, mj 


(3 ^ +4 ) G -1 r~j 

(2?,. sinh c .l)x(- 7-.- g .) j (EKNWN) 
J J ( A m +1) Aj ng , 


|(2,-X j (->l) j )x(^^~j) (0) j (EALPH) +|(2i j )x(G^) (2fi.)x(G^) 


-l^j 


(EBETA) 


(A. 18) 


•:\ m +2) 

m 


9 2 

1G 2 ) (0) 

m ( A + 


'm lG 2 ) 

( A +1) in 
m 


i*) (0) (2/ j (-l) j )x 


-- r^- T + rtj-y) (2 S 

A i I o A j 3 


( A +2) G 15 

sinh b' .l)x(^-r™pr.--^ + G* 3 ) 
J < ^ m +1 > --j m J 


(A +2) G 13 , - _ ( >, +2) G 18 . 

ar-T^ + (2> j sinh ' i 1)x( 

J v ' v m ‘-J y m • j 


( > +2) G 13 -r 

C. cosh .\l)x( r - ^ T .-^ + G ®) | (EKNWN) 
J J ^ ' x rn * c: j J 
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G 6 . 


, (O) ( 2 ( "1 ) 3 )x( - =.~r i " T ) ' (EALPH) + \(2 r . tanh6 . 1 ) x ( -G 1 ^ ) 

L J I l^j ; j I 3 x 3 mj 


g 13 : 


< 2 fj )x< - tanFfTl) ! (EBETA > 

3 


(A. 20) 


t* ( 3 '\ +4 ) a 

|(0) (0) (O) (_-_.iIL-._- G_) (O) (2 )x( 


( \ + 2 ) 6 °. 
m mj_ 


m 


G 7 . 


,15 


o 3 


( 2 . ) x ( ■ 


- . ( "S +2) G . 

+ ri^ry) (0) (2 ^ )x( 'n\ iI W*‘T I + G mj } (Q) 

- - / m *■ j 




+2 ) G 18 
m nji 

+ D* f, ■ 

m °j 

16, i 
+ G . ) ; 
mj ! 

(EKNWN) 

+ * (0) 

(2 )x(= 

G U 

_ . JL2 __ 

sinh (' .1 




I l ( ^ j ) j 


(A. 21) 


(iii) Equation on Interface (Bonded Region) : 

f, 3 1, , 1, , .,3, , 3 2, , „4, ,„6 , , aI o^ 

(a X ) (ax ) (ax_) (a x) (aX.) (X . )x(—_ 

m m m m m mj I (X .) 

L ° 3 


ai ( x ,- a ) 
o j 

r 3 1, , ..1, , _..3 


,8 


( A j ) x (“ i £ T~T~ ) (X mj )x(aJ o ( ^ j a) (X mj )x(aJ o ( j a) 1 (SKNWN) = 

ai ( u .a) 


ai ( ./ .a) 


l(a^) <<) (<) c»X> ( <j )x( ^77ir> ( A ? x <~i 2 rq J T 

o 1 3 o 3 

(x Ij )x(aJ 0 <> "j a,) (X mj )x<aJ o ( * j a)) <X raj )x< '’ x o (C - j a)) 


1-1 (>>\ J 3 ) 


1 , ( -\.- a ) 


(X?.)x(aY ( $ .a))‘j (EKNWN) + f(X^ . )x(- ~r^-r ) ( ^)x (•---■• /- ■■ )! 

nig o j ; o m 3 


3 ’ 1; 
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9 J-iCp.a) J-CC.a) . 

(EALPH) + j (X raj )x< 5sa iJ- T ) (EBETA) 

v r J Y j ^ 


6 1 ^ i ^ m 1 -t ( ^ . 3- ) j. 

(X mj >X( l7TTj'r ) < | (SALPH) 


"9 J i ( Fi a) 
— 1 


. 10> , J 1 ^ I" j a) ') 

(x mj )x( iiHrr.T ) ) (SBETA) 
i j 


, v A_ . j r — f — r 

mj cosh .a .1 
L J V j 

(A. 26) 


) 


4 ( /> j-+2) „ „ 8X 

L ( “ ri^Trr a x m + 3 Tt~tit ) (0) 


2, r 2 2 \f 4, 

7^+1 a ~ 3 TJT+TT V 


2 ( ' +2 ) 

zv ,. f +*> 2 3 2 1 

(„ ______ ax + 

( 4-1 ) m 


5 (2 A f +3) q -j_ 

3 ( '> _+lT X m + 4 ( 37+17 a X m } ( ' 


3 2 3 

2< A £ +l) a X m 


1 5 3 ( ^f +2) 4 1, , aI l ( <i a) 2 ( i f +2> I Q ( < i a) 

2( A f +n X m + 16 ( -A f +1) a X m ( I q ( <7) + r>“5TT I Q ( . ) 


m, ,„6 N/ aI i< x -j_ a ^ 2 ( 'A f +2) I n ^ *7 -i a ) 


x(»Z..-) (X-.)x( 


^ • > ' + ( Mirx/irn 


• T / '1 / 


J O 


<3j )x(aJ l< ■- j a) 


^'m* 2 2J o ( ' j a) 


m 


+ 1 * 


c 7 

mj 


(X„_. )x(aj, ( b .a) 


A m +2 2J ( S .a) 

P — — J ) ; (SKNWN) 


^ +2 ) O 9 Q 

I ( ^ rj2y2 8 

> ( 2 +1 ) X m + 3 ( ;s +TT 

m ' ' ra 


m 


2 2 


2 A 


2 <> m +2) _ 2 .. 3 


(0 > <- O^IT <- Ta^rr a x m + TA^Tr 


2/>t m +3 41. , ..«!<<£> , 2( ^ m +2) VA 81 , 
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6 % / aI l ( ^i a) 2( ' > ' m +2) 1 0 ( -<n a) 




I o ( ->: j ) ( A m +D 


J o 


^Y~) (X m ^)x(aJ 1 ( ^a) 


7^4-2 2J ( o .a) „ +2 2 J ( o .a-) 

'T+l- _) (X mj )x<aJ l ( '' j a) - — = L “> 


'’3 


m 


8 


\+2 2 Y (c .a) 


m 


(X^JxUY.U.a) 


° 2 ) (XL .)x(aY 1 ( < .a) 

mj ± j 


'\ +2 2 Y (>. .a) _ 

“T* ^ )1 (EKNWN) 

A m + <- j A 


r m 31 (-x, .-a) , 1(2 .a) 

• ( *=■ L. ^ ) x ( =~r-— ) X 6 .(-^ ' 

l. 3 I 1 (o< J mj I 1 (, K ) 


— )j 


in J ( n - a ) 

p ^ \ / O ' 1 


(EALPH ) + l(-X mJ )x( eoshp . .! 

¥ J 


> (EBETA > 


.'7 m N / I o ( i a \ 

?- A j )x< 5rr^r ) 


el I ( -V -a) 

(X*_ . )x(r~7 3L_) ! (SALPH) 


.10, 


mj 


'Ii < * j) 


J ( l" .a) 

° LL 


! ( “ X mj )x( cosh fi . 1 ' 


J ( f' : .a) , 

J— -) i (SBETA) 


( ” x mj )x( sinhp T 


(A. 27 ) 


< 2( Af+ 3) 2 1 8 ^f 3 ( 2 (>+l)X 1 ) < ( Af - 2 - 1 

'’F I X £ + l ) X m + ( » f +l) V 'F (2( A f +DX m ) A f ( /if 


A j; +1 X m 


+ 


2 ( 2 >\ f +l) 3 


( ) m 


xZ) 'U( 


2 ( A f +3) 2 2 8 A £ 


F v ( ;l+D 


ax + - 


xl) 


'm 3 ( ^ £ +l ) m 


S*< 


3 A 


2(2/1 f +l) 4 

a X„ + - 4 tt~ X ) ( 2 i' tX.xZ, ) C 


( ^ £ +l) * 3 ( L f +rr m' 


6 w aI l (c <j a) 


F ( " v j mj I o ( u .) 


A,- ( ,1^+1) *I Q ( X n . ) 


m w aI l ( < i a) . 
) (2 S y • ■£ • ) V~ — ; -H— + 


I CL .a) 

° 1 2 . 


J J i o («j) yj < ^ £ +Di 0 <Kj) 


4 - 
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7 1 J ( 7 .a) 


) (-2-^ % .X 8 .) 

>F j mj 


(aj i ( S .a) 
1 1 


( .x f+ D- . j 


J ( -l -a) n 
- g ~-r~ J — ) i 


(SKNWN) 


- 2 ( ;v +3) _ . 

■ ( _ S =2 1 

( V + 1 ) X m 
' m 


+ TS^+ry x m ) < 2 <X +1 > X m> <• 


a V + x 3, 

■' rr +1 m ( ,V m +1 > m 


( 


^ m +3) 2 2 m 4 , 6 w ^lj 

r ; r+TT X m 3( A m +l ) V ^j X mj )( -‘f o (7 y ; j 


aI l ( j< n - a ) 


77 r + 


I ( X ..a) 


ali U' .a) 


1 V'r', . 1 I o ( <i a) 1 

<Y m +u j-vTr* 


J ( A .a) 


( " 2( j X mj )( ' aJ l ( : j a) ( ^ m +l)*“ £ 77 X “ ) ( “ 2 C J X mj )(aJ l ( ^ j a) 


rrw-if'"'’ { ~ 2S i x Ij )(aY i <c j a) " tjTit 


Y ( £ .a) 
,_2__LL__) 


8 


2 kj x m j )(aY l ( ^ j a) 


( A +1 ) * 
m 


Y < 7 a) 

~ 2 -— — 1 ) 1 (EKNWN) + 


I (x .a) 
]( 2 V X° . ) -2 — t—7 — r 

t m j ; * ^7777 


i Ii(<.a) I (-X .a) 


7 ^* 1 ^ 7777 “' v ^j kr^TTTT 


I, ( 7 .a) 


7p:-5rr^ 1) ’! (EALPH) 7. <X "J :V cosh Pj 1)(J ° (jS j a) - 77 al ) 


in ~ 1 ^ i a ^ "] 

(2x mj *- /smh r' j D(J 0 ( p ja) - 


(EBETA) 
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, I o ( -< i a) I x (x .a) 


p N F ^j x mj )( i' 1 rpr “ ^Tiip"x j ) ) 


(2 


I .a) 

, , nu . c>< j 

" r ^ 'J B'-Ki* 


x ^ ^ 


. 2 ) j (SALPH) ~ i (2 .? ./cosh V.1)(J ( .a)- 

ai i o( . ) j , i j / j o i j 


(2 ^./sinh f.l)(J o ( p.a) -4-^ 


r j' 3 ' " o’ ] ] a : 

) yep I (SBETA) 


^Ap<x 9 . 

F mj 


*' j 


1< 


JO 


mj 


(A. 29) 


8(;\,+2) „ 

V" rxr+TT "V <°> 


2'A 


T A^r+l) m 


ax“) P(- 


3 A £ +4 


3 1 

Ac+1 a X m “ 


4( ' Xf+2) aY 3 ^ -< ( 3( ^f +2) 3 1 + 

TP .+1) " aX m 4 ( \.+l) a X m + 


3-v 

( >. f +l) aX m ( “ 2 ‘f X j A j )x 


( 


jp • *** / **«■ * V ^ f 

ai ( ,>„/ .a) ( A + 2 ) Ii(‘X -; a ) 

+ < > f +D j) > (2 ^f ^j )( l 0 (^j) 


, al ( j .a) 
6 w o a j_ 


,('A +2) I-i ( T .a) 

Trprr^- i7^r > ' <2< f + 


8 w __ *f +2 J l<£ i a > 


A,+1‘ 


O 


) 


(2 t _ S -X 4 ) (aj ( .a) + — — 

■ F °j m] o j ,X f +1 


p + 2 J.U -a). 


) 


(Sknwn) 


r, 8( ' X m +2) _.,2 


2 A 

•’ m 


u-r^rr^ (0) ( T3^IT<> <T^ 


3 ^ m 44 3 1 


a X_ 


m 


4( ^m +2) „„3, , „ . m., aI o <j( 1 a> * , 

Tyrr*. 1 ( - 2 =<j *j )( TX*Tr + TTjzr~-ijr^ 


( V- +2) I, (x.a) 
m 1 }) 


m 


J o "-J 


al„( K,a) (> m +2) I 1 (^.a) 


(Oni v b ) ( ° - -TL- + '_g .J- IX, 

U< i X ™j n I 0 (^j) (A ra +1) < J ' I o<-<j> 


p~) ( 2 ; X® . ) ( aj ( 2 r .a) 

3 mj o j 
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A m +2 J l ( t-j a ) 7 _ 'X m +2 J i( <: M a) 

Arn j J ra J ° 3 ; j ' 


8 *- Y- ( c a) 

(2 e j x mj )(aY o ( f ja) + TTXrr •— --- J > < 2 <■' jVj > (aY o ( ^ j a > + 


^n, +2 Y l (; i a >X 

T3' — rr — > 


(EKNWN) 4* ! (-2* . ,<^)( r *L_L ) 

» J J •£ i \ >> . J 

»- - 1 J 


, „6 *1< ■<.,*>. , 


<2 *f£j )( str^f > J <EALPH> + | <2X m? f/coshhjlHj^ f,.a>) 


( 2X 


9 ,7 


m. ^l^'V - a ^ 


mj 


(,'j/sinh j J l) (J x ( Pja)) j (EBETA) =!(»2- SF ^ j ) 


(2 5F-’'.j X m3 )( l7r3^r ) j <SALPH) -4 <2 <F X mj f/ cosh 


(2 ^F X mj f'j/ sinh pjl)(J 1 ( | (SBETA) 


(A. 28) 


( i v ) Equations on interface (Debonded R egion ) 
f 2(A f+ 3) 2 8 

[ ( r>^rr a + 


f~ r '" _ 2 w £ 2 


(T^lT z rn } (2( -'f + l ) (- “f~T=T a 2 


W+T ° + 

£ 


2 ( 2')\ £ +l ) 2 

( 3\ r: +1 ) m ^ 


3,2 


2 ( ^f 4-3 ) 2 8/X f _3, , “ 2- f 2 

( >,+1) a Z m + 3( Ar+1) Z m ( ~ “TL+T a z m 


2 ( 2 /\ £+ 1 ) 3 


, aI l ( A i a ) I Q ( X.-a) 

H . mnl L. u _L» — I , mi ITT 


+ 3Tli^4TT" z m ) (2 y 'j COS ^j Z m ) (“l^Tx~T " T "V7T l f +l)l Q (~ j ) 


) 


( 2 o( . sin ,v .z ) ( 


al n (x ,a) I ( x ,*) 


'! v i 


2 Li 


•jm ,v I o (^ j ) T o< j I o ( v .)( A f +1) 
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(-2 fcj cosh b j z m ) (aJ x ( Lja) - 


1 J J ^ 

1 °_™JL__) 


J (F, ..a) 


(=■2 hj sinh j 2 m )(aJ 1 (f, ^a) - 1 — ) f (sknwn) 


+l<2 -< s cos f-h.j y> <2 S sln 


i ( y - a ) 
1 


1,(2^) 


cosh p z 


x^.xp I (SALPH) +| < 2 ? i cbihjjrxv^ j a) 


J i< Pj a ) 


sinh ^ .z 


i j ; j 


j i^ r i 51 -* ■> 

-r— r 3 ) : (SBETA) = (^) 

■ a i 


f j 


(A. 36) 


r 8(A f +2) 


2 'A 


rxprr a V <°> Acrr a V <- 


3 Af+4 3 4 ( f + 1 ) 2 


-+1 


T - or-,- az ) 

( 4 f +l) m 


3 ( ’A f +2) 3 


3 A 




aI °in a) 


■4T17P1T a ~ + rTT^TT az m> ( - 2 A sln -'j z m >( -Trr7tr + 


o 'S j 


( / £+2 ) ( •«( ja ) 

TTpiTTj-i77j) 5 


ai ( y .a) 

( 2 x . COS y‘ . Z ) (——§■—•; — + 

3 3 m I Q { 'X j ) 


( ^ £ +2) I x ( \ j a ) 

T X £ +l) rXj ^ *■ j ^ 


(2 A . sinh ' .z )(aj ( p .a) + . 

3 3 m o j A ,+l 


:+2 J, ( V. .a) 


1111 


■ 3 


(2 F cosh l. z )(aj (f a) + 
J j m o j 


;i f +2 J 1 ('. ^a) 


/\ £ +1 c j 


) (SKNWN) + 


r; 


I 1 C< ,a) 


I-, ( < .-a) 


sta •<i 8 m )< itr^r ) <2< j 008 ft >j 


(SALPH) 
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sinh f _.z 


+ (2 


Pj ssr^r )<J i ( pj a) 


cosh U ,z _ 

(2 Pj sinh £ “!i m) (J 1 ( i 1, j a) * j 

lh j J j 


(SBETA) = (O) 


(A. 38 ) 


20 + 3 ) „ 8 . 

m _2 m2 


( > m +T ' 


L ' m 


2 ( / „+ 3 ) 




: m 


a + 


a 2 z 


'A +1 m 

s m 


O (2( A +D ( 


M m 


2(2A +1) 

^ m 


m 


\ -hi 
^ m 


a + 




m 


- 3 — 2 2 ) 

+1 m 


8 A 


ry m -?rr a z n + z m> ( 2 -<j cos <- i-jrfr + 


al ( ,,< .a) 

-Ju_™OL 


' m 


1 ( < - a) ai a) I ( >( .a) 

+i)i. ( .< 7J ) ( 2 :< i sln j 2 ™’ ( i c ( ^ j) + .~J(; lm +i)ioKj) 1 


,X j (A m + 1 ' I o ( - < j ) 


(-2 C } cosho ,z m )( a j 1 ( t j3 ) - T f ~ hTiO (_2S 1 slnh ^ j z ra )x 


J _c< ; ^ a > 

M m • i-j 


(aJ i ( ‘ j a) - — ( " 2i ’i cosh -i z m )(aY i ( >■ j a > 

' / ‘j >yi / 1 -*■ “i •*' 


rfO-f o ( - 2 ■' j sinh ? 'j 3 m )<aY i< i- j a > - 


Y-( £ ,a) 
nO - j 


p. I ( to \ . a ) I 1 ( ^ .a) 

ceknwn) + : < 2 * cos -x z ) - — iO'-ry) 


* ^ ,a) I. ( - ,a) 


(2 U . sin .Z ) (-— -ly- - -0 

J 3 rn !-!_( :*■ j ) «X k jai 1 ( ^ 


i: i 


) I (EALPH) + 


|(2 fljtcosh fj z m / COSh f j 1 )XJ o ( Pj a) ' 

) ' J 


sinh [f 4 .z 

< 2 Pj a 5 5K - ^r )(J o ( f 3 a) 


J x ( ,- a ) -t 


Pj 3 


^ ); (EBETA) = (O) 


(A. 35 ) 
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8 (,'\ ■ + 2 ) 2 h 3 > +4 , 4 ( /. + 2 ) 

(- 7- — — az ) (°) (— -2— az ) (_~ 5 L_ a 3 - 7 — — — az 2 ) 

( m >„+l m' >. +1 ( /\ +1 ) m J 

' ' m m 


■ m 


* m 


( “2 x . sin ^ .z w ) ( 


aI ° ( K i a) j. /N m +2 J l ( < j a ) 


m_ . j ( 7 .y 

j O - ' j 


+ 


m 


+1* < ,1 (x .) 
m J o ^ j 


ai^( ' ,a) A m +2 I x ( ; . a) 



'V + 2 j l ( ^'i a) 

( 2 .: sinh : 3 z m )(aJ ( > a) + -^.—-2—) 

j j • m ' j 

A m + 2 J i ( A .a) 

(2 cosh > j 2 m )(aJ o (?: a) + 

j j j m • j 


^ Y- ( c .a ) 

(2 sinh J) z ) (aY ( .a) + — I T .— — ^ — ) 

j J ^ m j 


(2 6 . cosh ; .z )(aY ( ; .a) + — — — — 
J 3 m o 3 y \ m +l \ 


\ +2 Y 1 ( c .a) 

“ — - * -L- -) | (EKNWN) + 


c 3 


I 1 ( 'A . a) 


k-2^ sin ■<;,»„> (i^ r 


-) (2 t cos o' .z ) ( 


^ I ( X ^ a ) 1 

^ ) j 


. wo c/ . * > \ - 7 ■ S ; j 

J J m j 


sinh f .z 

(EALPH) + j (2 V J PjS)> 


cosh if- .z 


< 2 rj id 5 T^r )(J i ( ?j a)) ! (EBETA) = (0) 


(A. 37 ) 
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APPENDIX III: FLOW CHART 


Evaluate Fourier and 
Fourier =Bessel 
coefficients 


Bessel's function routine 
-Trapezoidal integration 
routine , Block Data 













STIC MODEL 






Fig . 3.2 Variation of nondimensional sheaf stress 
with nondimensional length, 



Fig. 3.3 Variation Of nondimensfbnat rgcha 1 displacement 
■ ''' with nondthnensional length.-', 



